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The relaxation rates to the invariant density in the chaotic phase space component of the kicked 
rotor (standard map) are calculated analytically for a large stochasticity parameter, K. These 
rates are the logarithms of the poles of the matrix elements of the resolvent, R{z) — {z — U)~^ , 
of the classical evolution operator U. The resolvent poles are located inside the unit circle. For 
hyperbolic systems this is a rigorous result, but little is known about mixed systems such as the 
kicked rotor. In this work, the leading relaxation rates of the kicked rotor are calculated in presence 
p^ , of noise, to the leading order in 1/y/K. Then the limit of vanishing noise is taken and the relaxation 

rates are found to be finite, corresponding to poles lying inside the unit circle. It is found that 
the slow relaxation rates, in essence, correspond to diffusion modes in the momentum direction. 
Faster relaxation modes intermix the motion in the momentum and the angle space. The slowest 
relaxation rate of distributions in the angle space is calculated analytically by studying the dynamics 
of inhomogeneities projected down to this space. The analytical results are verified by numerical 
simulations. 

PACS number(s): 05.45.-a, 05.45.Ac 
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O' For chaotic systems specific trajectories are extremely complicated and look random [|l|. Therefore 

_^^ . it is natural to explore the statistical properties of such systems. For this purpose the evolution 

Qs^ . of probability densities of trajectories in phase space is studied |@-^. For chaotic systems the 

"^ I probability densities, approach an equilibrium density that depends only on the system and not 

>^- on the initial density. For hyperbolic systems (A systems), like the baker map, the relaxation is 

exponential. For such systems the existence of the relaxation rates was rigorously established and the 

O [ relaxation rates are the Ruelle resonances 1^-01 • To study these rates it is instructive to introduce 

(--| ■ the evolution operator of densities that is sometimes called the Frobenius-Perron (FP) operator. 

O ■ Relaxation to the equilibrium density is studied traditionally in statistical mechanics. In particular 

>• . for particles performing a random walk in a finite box, relaxation to the equilibrium uniform density 

k> [ takes place and is governed by the rate related to the lowest nontrivial mode of diffusion equation. 

^ It is known that for the classical kicked rotor, described by the standard map, diffusive spread in 

phase space takes place for a sufficiently large stochasticity parameter PJ^. Therefore it is natural 

to study the Frobenius-Perron (FP) operator for the kicked rotor and to compare it to the diffusion 

operator, a comparison that enables to study some aspects of chaotic dynamics in the framework of 

statistical mechanics 0]. The kicked rotor model is a paradigm for the chaotic behavior of systems 



I. INTRODUCTION 



where one variable is unbounded in the phase space. For such classical systems diffusive spreading 
takes place. For their quantum counterparts it is suppressed by interference effects, leading to a 
mechanism that is similar to Anderson localization in disordered solids [1TT|-[1^ . 
The kicked rotor is defined by the Hamiltonian (in appropriate units) 

n = -.P + KcoseY,Sit-n), (1) 

^ n 

where J is the angular momentum, 6 is the conjugate angle {0 < 6 < 2tt) and K is the stochasticity 
parameter. Since the angular momentum between the kicks is conserved, the equation of motion 
generated by (|l|) reduces to a map, known as the standard map 



6 = 6 + J (2) 

J = J - Ksin6, (3) 

where 6 and J are the angle and the angular momentum before the kick, while 6 and J are these 
quantities just before the next kick. For K > Kc ^ 0.9716 diffusion in phase space is found, and for 
large K the diffusion coefficient is given by an expansion in -^ as: fl^ 



DiK) = ^{l-2MK) + ...). (4) 

To be precise it was shown that after a large number of kicks, n, the variance of the momentum 
behaves as: 

((J-(J))2)~2Dn, (5) 

where (■ ■ ■) denotes an averaging over the angle initial distribution, and D is given by (^). 

It is assumed that the system evolves in presence of finite noise and the limit of the vanishing noise 
is taken in the end of the calculation. The noise is required here in order to get well defined results. 
It leads to escape from the accelerator modes and other stable islands. Accelerator modes, where 
the angular momentum J grows linearly with time, are found for values of K and the initial values 
{6o,Jo) of the angle and the angular momentum, that satisfy Ksin6o = 2ttIo and J = 2ttI where / 
and /o are integers. In such a situation at each step, J grows by 27r/o, as is obvious from (|) and (||), 
namely its growth is linear in time. For some values of K the point {6q,Jq) is stable and also for initial 
conditions in its vicinity the momentum grows linearly. This differs from diffusion, that takes place 
in the chaotic component of phase space. For trajectories in the chaotic component of phase space 
noise avoids long time sticking in the vicinity of islands of stability [15]. In numerical calculations 
without noise diffusion is found for K > K^ for trajectories in the extended chaotic component for 
large values of K, but also some exceptions were reported JlSl. The diffusion coefficient @) was 



calculated in presence of finite noise (in the long time limit) and the limit of the vanishing noise 



can be taken in the end of the calculation [Q. It describes the typical spreading of trajectories in 
the chaotic component. Since the kicked rotor is a mixed system, as is the case for most physical 
examples, the rigorous mathematical theory for relaxation 0,|^,0] does not apply and one has to 
resort to heuristic methods. 



In the present paper [T^] the Frobenius-Perron operator will be calculated for the kicked rotor on 
the torus: 

(0 < J < 2txs) (6) 

(0 < < 27r), 

where s is integer. This is reasonable since the map (0-§|) is 27r periodic both in J and in 6. The 
operator is defined in the space spanned by the Fourier basis 

11 / kj\ 

<Pkm = {J6\km) = — —i= exp{im6) exp i — . (7) 
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Note that the functions (pko form the basis of eigenstates of the diffusion operator in the angular 
momentum J. The FP operator for an area preserving invertible map, 

X = M (x) 



IS 

Up{e,j) = p{M~\e,j)). (8) 

It was studied rigorously for the hyperbolic systems and many of its properties are known 
PJ3|j5|,p!^,p!^ . It is a unitary operator in >C^, the Hilbert space of square integrable functions. There- 
fore its resolvent 

1 1 °° . 

R{z) = ^ = - E U-Z-- (9) 

z-U 2: „=o 

is singular on the unit circle in the complex z plane. The matrix elements of R are discontinuous there 
and one finds a jump between two Riemann sheets. This results from the fact that the spectrum is 
continuous and infinitely degenerate 0. The sum (|^) is convergent for \z\ > 1, therefore it identifies 
the physical sheet, as the one connected with the region \z\ > 1. (This is analogous to the sign of 
the small imaginary increment in the energy that is used in the definition of the Green function). 
The Ruelle resonances are the poles of the matrix elements of the resolvent, on the Riemann sheet, 
extrapolated from 1^1 > 1 [0. These describe the decay of smooth probability distribution functions 
to the invariant density in a coarse grained form 0,0. Even a smooth initial distribution will develop 
complicated patterns as a result of the evolution of a chaotic map. The Ruelle resonances describe 
the decay of its coarse grained form to the invariant density. In spite of the solid mathematical theory 
there are very few examples where the Ruelle resonances were calculated for specific systems. They 
were calculated analytically for the baker map where the basis of Legendre polynomials was used 
T7| and various of its variants 0. The Ruelle resonances were calculated also for the "cat" map 



and some of its variants [^]. Blum and Agam applied a variational method for the calculation of 
the leading Ruelle resonances of the "perturbed cat" map, and the results were verified numerically 
||T9| . In addition they calculated the leading resonances of the standard map with s = 1 for various 
values of the stochasticity parameter K. The leading Ruelle resonances for the kicked top were 
calculated by Weber, Haake and Seba |2^ with the help of a combination of a cycle expansion 



and numerical calculations. The resonances mentioned above are not related to the spectrum of the 
Liouville operator that is confined to the unit circle because of its unitarity. 

In the present work, the FP operator is calculated for the kicked rotor. Here the classical evolution 
operator, for one time step, can be written in the form 

Uj^^ = 5{e-e-J)5{J -J + KsinO), (10) 

and its operation on a phase space density p is 

Up{e,j) = p{e-j,j + Ksm{e-j)). (ii) 

To make the calculation well defined noise is added to the system. It is shown that the noise acts 
effectively as coarse graining and the resulting FP operator is not unitary (see also [^). For large 
stochasticity parameter K, we show that the slowest relaxation modes in the limit of infinitesimal 
noise are the modes of the diffusion operator in the angular momentum space. Also calculated is 
the slowest rate of relaxation in the angle space. The approximate analytical results are tested 
numerically. 

It is understood that the noise is kept finite when the limits of large K and s are taken and then 
the limit of zero noise is taken. The natural question is whether it is possible that this description, 
which was established only for hyperbolic systems, holds also for mixed systems. Clearly, for mixed 
systems it can be only approximate. It holds for large values of the stochasticity parameter K since 



then most of the phase space is covered by the chaotic component. For smaller values of K the 
weight of the regular regions increases. In such a situation, in the limit of increasing resolution the 
resonances related to the regular component are expected to move to the unit circle in the complex z 
plane, corresponding to the quasi-periodic motion, while the resonances associated with the chaotic 



component stay inside the unit circle. This was found by Weber, Haake and Seba ||20[ for the kicked 
top that is a mixed system. 

How is the FP operator related to the quantum mechanical evolution operator? It was shown 
numerically for the baker map that if both operators are calculated with finite resolution they 



exhibit the same Ruelle resonances PTl. In this calculation it was assumed that the phase space 



coarse graining tends to zero in the semi-classical limit ^ — ;> 0. It was shown by Zirnbauer |^ 
that some noise is required for a meaningful definition of the field theories introduced to study level 
statistics for chaotic systems P3[. This noise affects only quantum properties, therefore the resulting 



ensemble has the same classical FP operator. The localization length of the kicked rotor calculated 
from this field theory [^] is related to the classical FP operator. This operator is analyzed in the 



present work, clarifying some issues of that work. The results hold only for typical quantum systems, 
since the noise introduced in the present work as well as the noise required for the stabilization of the 
field theory p2| washes out the sensitive quantum details such as the number theoretical properties 



of the effective Planck constant [ Pq] . 

The Frobenius-Perron operator in the basis (|^ in the presence of noise is defined and calculated 
in Sec. 2, its Ruelle resonances are obtained within some approximations in Sec. 3 and their regime 
of validity is tested numerically in Sec. 4. The results are summarized and discussed in Sec. 5. 

II. THE EVOLUTION OPERATOR OF PHASE SPACE DISTRIBUTIONS 

In this section the evolution operator of phase space densities of the kicked rotor in presence of 
some type of noise is derived. The noise is added to the free motion part (^. In the absence of noise 
the phase space evolution of a distribution / is given by Liouville equation 

df _ df ,df -df _ 



If noise that conserves J, and leads to diffusion in 9, is added to the free motion, equation (|l^) should 
be replaced by 



91 jdl_ a^^f 
dt 89 2 09 



^ + ^^-^^ = 0. (13) 



where J = 9 was used. It can be written as: 

d 
~dt 

where the operator A is: 



% = Af, (14) 



The complete set of its eigenf unctions is given by (fm = ^= exp(i?Ti6'), where m is integer. The 
operator we need is Unoise = e^, and explicitly 

{J'9'\Ur,o^se\J9) =Y.y^U0')^*m{G)exp{aJ6{J - J'), (16) 



where the «„ are the eigenvalues of the operator A, namely A\iprn) = c^ml'^m)- Obviously 






leading to 



Ctm = —imj —771 (17) 



{J'0'\Uno^se\Je) = E ^ ^Xp (77n{e' - 9 - J) ~ ym^^j 5{J - J') . 



The (5-function in momentum reflects the fact that the noise does not affect the momentum. 
The matrix elements {k27n2\U\kimi) of evolution operator U in the Fourier basis (|^) will be 
calculated in two steps, first the contribution of the kick, and then the one of the free mo- 
tion with noise will be calculated. According to (j^) and (p!TD, the kick transforms the state 
(6*, J\ki7ni) = -^—^ exp(zmi6') exp (i^j to the state 

exp(zmi6') exp i — (J + Ksi7i9) = {9, J\UK\ki7ni). (19) 



V27r V27rs V ■^ / 

Adding the effect of noise yields the matrix element in the mixed representation 

r2n r27rs 

{J9\U\ki7ni) = d9' dJ'{J9\Unoise\J'd'){J'd'\UK\ki7ni). (20) 

JO Jo 

Its transformation to the basis (0) is calculated in App. A and the result is 

{k2m2\U\ki7ni) = Jm2~mi f -^— J exp ( — -ml J Sk^-kumis- (21) 

For (7 = 0, using ( ^I]) one can verify by a straightforward summation that WU = I, therefore the 
operator is unitary as required. 

Some of the eigenfunctions of U in the limit a = are easily found. It is convenient to use the 
representation (p^Of ) of U. We guess an eigenfunction of the form 

F{9, J) = 6{9 - 9q) E exp(zgj/s)(5( J - 27r/) (22) 

with q integer satisfying 1 < q < s. These are functions localized on accelerator modes representing 
linear growth of the angular momentum with time. To check that these are indeed eigenfunctions 
we note that 

UkrF{9, J) = 6{9-J- 9o) E '^(^ + Ksi7i9o - 2-kI) exp [i- ( J + iTsin^o)) (23) 

I \ s J 

taking 9q so that 

Ksin9o = 2nlo, 

where /q is integer yields 

UKRF{9,J) = e'^-''''^/'F{9,J), (24) 

since the RHS of (p3D does not vanish only for J = 2nl. The eigenvalues e*^'^'^'"/* lie on the unit circle 
and become dense as s — > cx). There are more eigenfunctions of this form located on other periodic 
orbits |26l . 



III. IDENTIFICATION OF THE RUELLE RESONANCES 

The purpose of this section is to calculate the Ruelle resonances for the kicked rotor with the 



help of Frobineus- Perron operator (0). The calculation will be done for finite noise a and then the 
limit o" ^ will be taken. The Ruelle resonances are the poles of matrix elements of the resolvent 
operator i? of (H), 



Ri2 = {kimi\R{z)\k2m2) = kiini 



z-U 



k2ni2 , (25) 



when analytically continued from outside of the unit circle in the complex plane. It is useful to 
introduce the operator 



that is related to the resolvent by 



and 



r^{z) = ^, (26) 

1 — zU 



^-E! (-] = R{z) (27) 



z \z 



^-R {-] = R'{z). (28) 



z \z 



The matrix elements of R and R' satisfy similar relations. Continuing the matrix elements of R{z) 
from outside to inside of unit circle is equivalent to continuing the matrix elements of R'{z) from 
inside to outside of unit circle. The last continuation is easier to study since the expansion 

R'{z) = - ^ = E^"^" (29) 

I- zU n=0 

is convergent inside the unit circle, because 

\\zU\\<l. (30) 

The resulting matrix elements are 

oo 

R',2 = ikim,\R'{z)\k2m2) = ^ (^nz"", (31) 

n=0 

where 

a„ = {kimi\U'^\k2m2). (32) 

Through ( ^Tf ) and ( ^Sf ) this expansion is related to matrix elements outside of the unit circle. If z^ 
is a singularity of Ru then 1/zc is a singular point of R^2- Consequently the first singularity of the 
analytic continuation of Ri2{z) from inside to outside the unit circle gives the first singularity one 
encounters when analytically continuing Ri2{z) from outside to inside the unit circle, i.e. it is just 
the leading nontrivial resonance. This is the most interesting resonance determining the relaxation 
to the invariant density. The first singularity in the extrapolation of the matrix elements of R! from 
inside to outside the unit circle is determined from the fact that it is the radius of convergence of this 



series. Moreover according to the Cauchy-Hadamard theorem (see [^) the inverse of the radius of 
convergence is given by 



r ^ = hm sup {/knl- (33) 

n^oo i V ' ' 

If a„ ~ cjr"^ we may say that the radius of convergence is the asymptotic value of an-i/cin- This 
is the basis for the ratio method for determining the radius of convergence. The resonance that is 
closest to the unit circle can be identified from the radius of convergence. 

We turn now to calculate the coefficients a^. First the matrix elements {kQ\E! {z)\k{]) will be 
calculated. For these the expansion coefficients are: 

a„ = (fcO|f/"|A;0). (34) 

Introducing the resolution of the identity, 

«n=E E- E mU\kimi){hm^\U\k2m2)...{K-im^^i\U\kQ), (35) 

fci,nii fc2,m2 fc„_i,m„_i 



and substitution of ( plj) leads to: 



■n — ^ ^ ■■■ 2^ -'0-mi I j Ok-ki,0-Jmi-m2 \ I GXp I — — m^^ I Ofcj_jfc2,mis l<JDj 

ki,mi k2,m2 k„—i,m„-i \ / \ / \ / 

(k,K\ (y ,\ (kK\ (y , \ 

'Jm2-m3 I I exp I 1712 I '^k2-k3,m2S---'Jm„-i~0 I I 6Xp I ^^n-1 I '^k„-i~k,m„-is- 



Summation over the ki yields, 



EE-E^o-™.f^) (37) 

■mi m2 m„-i \ ^ / 



Jmi-ma "^1-f^ exp — m^ Jra2-m3 ("^1 + Vfl2)K CXp — TTlj 



■Jm„.2-m„-i (mi + m2 + ... + m„_2)A exp — z-m„_2 



'kK\ I a- 



"^mn-i-O I — ) exp ( —fTT'n-l ^mi+m2+...+m„^2+m„-ifi- 



Thus in order to obtain the expansion coefficient a^ we should perform summation in ( P7D over all 
integers subject to the constraint 

mi + 1712 + ■■■ + ^TT'n-l = 0. (38) 

We are interested in the limit of large s and K. The limits are taken in the order 

(1) s -^ oo, (2) K -^ oo, (3) a^O. (39) 



Finite a is required to assure the absolute convergence of the series. Therefore (pTl) is summed for 
finite a and the hmit a — *^ should be taken in the end of the calculation. Having this limit in 
mind the leading term in K/ s and 1/\/K will be identified. It will be assumed that the mode that 
is calculated is sufficiently low so that 

< kK/s « 1. (40) 

Each term in ( P?] ) is defined by the string 

(mi,m2,...,mi,...,m„_i). 

The leading contribution 

results from the string where all ttij vanish. A nonvanishing rrtj results in a Bessel function with 
a large argument, since K/ s « 1 and K » 1, and therefore it leads to a factor 1/yK in the 
contribution to a„. Let mi be the first nonvanishing ttij and mf the last nonvanishing one. The first 
factor in (0) that is not Jo(^) is J-m, (^) and the last factor that differs from Jo(^) is 



Jrrif I (mi + m2 + ... + mf)K \ = J^\ I . 



(42) 



Since J„(x) ~ 2^ ^^ small x and J-n{,x) = (— l)"J„(x) the contribution of the terms between i 
and / is of the order 

C (^) , ,43) 

where C is the contribution of the factors with ttlj that are not the first and last ones. The first 
factor after mi is 

Jm,-m,+-, I (mi + ... + mi)K I = Jmi-mi+1 I rTliK I 

and the last factor before mj is 

( kK \ ( kK \ 

Jmf.i-nif ( (mi + ...+mf_i)K\ = Jruf^i-nif ( {mi + ... +mf_i)K\ . 

The terms in between are of the form Jmj_i-mj (^ — ^j^ji where Mj = rrii + rrii+i + ...rrij 7^ 

that are of the order \j\fK. Therefore the largest contribution from a string m,, m^+i, ..., mj is from 
the shortest string, namely / = z + 1. Because of (|38D m^ = — m,j and because of (^) the leading 
contribution is from the string mi = —mif = ±1. The resulting contribution is 

C = J^^.m, ('^ - m,K\ = J±2 ('^ - {±1)k] ^ J2{K). (44) 

The string can start at n — 2 places, therefore the leading correction to a^^ is: 



.." ^ 2(„ - 2Ur (^) J. (^ - A-) Jl (^) e-" 



that is approximated as, 

a^y ^2(n-2) | 1 - 



eK^V^ 



4s2 



..(.,(f)e-. 



(45) 



The sum of the contributions (|4ll) and (H5|) is 



,(0) + ^(1) _ 1 _ 



2l^2\" / k'^K'^\'^~^ 



k'K 

4s2 



+ 2n 1 - 



4s2 



^^(-)f -1^ 



1 






(46) 



In the leading order — ^^5- ^ 1 and 
1 — . ■■> 



Therefore in the leading order 






1 



n-2 
lim = 1. 

n — ^00 Yi 



kK\\^^^ 



4s2 



+ 2J,(/<-) - e 



2 1^2 



Fir 



l-^(l-2J.(ir)e-^^ 



(47) 



The resonance closest to the unit circle, z^ is identified from ( |55D as the inverse of the radius of 
convergence 



Zk = ^ 



+ -r^'^J2{K)e-'' 



4^2 4^2 



or within this order of the calculation as 



2 7^-2 



1^'K 



Zk = exp 



4s2 



■{1 - 2J2{K)e-'' 



These are the eigenvalues of the diffusion operator with the diffusion coefficient 

D{K) = ^{l-2MK)e-'), 



(48) 



(49) 



(50) 



in agreement with the earlier results [14 



The analysis of the off-diagonal matrix elements 

a„ = {km\U''\k'm') 
is similar. We assume again ( |iD| ) and K » 1. Analogously to ( P7| ) one obtains 



kK 



J2J2--- J2 Jm-rm mK exp 



0" 



-m 



(51) 



(52) 



mi m2 mn-i 



J, 



'kK 



a 



mi— m,2 



(m + mi)K exp m^ J, 



'/tir 



— (m + mi + m2)K exp 



O" 



-rrir, 



• ••^m„_2-m„_i \ — - [m + mi + m2 + ... + mn-2)K \ exp I -ym„_2 I 

•Jm„-i-m' I I 6Xp I '^n-l I "{m+mi+m2+...+r?i„_2+m„_i)s,fe-fc' 

Because of the last 5— function the On 7^ only if (A; — k')/ s = g is integer. The leading contribution 
results from the string mi = — m, m„_i = q and all other m,j vanish. It is therefore of the form 

af = B Jr' (^) , (53) 



where 



'kK\ ^ (kK\ ^ fk'K' 



a 



2 



B = J2m{mK)J^rn ( '^^^^ 1 J-q ( '^^^^ 1 Jq-m' ( '^^^^ 1 exp(-(T^m^) exp ( — -q^ \ (54) 

that behaves as a^*^-* of (^Tj) in the large n limit. The leading correction is found from neighboring 
pairs wii = — ttij+i = ±1 as in the case studied before with a result similar to a^^-* of (^) in the large 
n limit. Therefore no new resonances are found from the off-diagonal terms with fc 7^ 0, in the order 
of approximation that was used. 

For s 3> 1 the diffusion modes in momentum space constitute the slow degrees of freedom of the 
system. However, the faster relaxation modes (or, alternatively the modes of a system with s ~ 1) 
cease to be angle independent. To evaluate the magnitude of such a fast relaxation rate within our 
perturbation scheme, we have to calculate matrix elements associated with relaxation of disturbances 
from the invariant density that involve functions from the angular subspace { \0,m,) } with m 7^ 0. 
Consider, therefore, the matrix element 

an = (Om|t/"|A;m'). (55) 

The equation corresponding to ( p6D is: 



^n= J2 J2 ■■■ J2 {0m\U\kimi){kimi\U\k2m2)...{kn-^imn-i\U\km') = (56) 

fcl,mi k2,m2 fc„_i,m„_i 

2^ 2^ ■■■ 2^ Jm-rm\ I exp I —m \ d-ki,msJmi-m2 \ I eXp I —mi \ dki-k2,mis 

ki,mi k2,m2 k„—i,m„—i \ / \ / \ / \ / 



— I exp I -y"^2 I 0k2-k3,m2S...Jm^^^-m' I — I CXp I " y"" 



Jm2-mz I exp I 777,2 I '^fc2-fc3,m2S---'^m„-i-m' I I exp I "^n-l I ^kn-i~k,m„-is 



and summation over the ki yields a nonvanishing result only if k/s = q that is an integer. In this 



case: 



mi 1712 ni„_i \ ^ / \ / 



Jrn2-m3((-"^-"^i-"^2)-f^)exp ( — ^"^2] ..■Jm„^2-m„-i{{-m-mi-m2-...-mn-2)K)exp ( — —ml^2 



10 



a 



2 



Jm„_i-m'(?^)exp I -ym^^J 6 



-m—mi—m.2 — ...—m„-2—m„-i,q- 



The result is independent of s. This is a sum over rrii constrained by 

m + mi + 1712 + ... + "^n-2 + mn-i = -g. (58) 

In every particular term in this multiple series, generally, we will have multiples of terms J^{MK). 
If M = and z^ 7^ 0, such a term vanishes while if both M and u do not vanish J^{MK) ~ -4^. The 
leading contribution is from sequences with the maximal number of factors Jo(0) = 1. To identify 
these, we denote Jo(0) = 1 by "1" and other factors by "x". In this way, to every term in ( p7D 
corresponds the sequence of n symbols: 

x*x*l*x*x*x*l*x*l*x*x*...*l*x*x. (59) 

A crucial restriction is that if m 7^ 0, two "1" symbols cannot be nearest neighbors as is shown in 



what follows. If TTi 7^ the sequence starts with "a;" as is clear from (|57|) . Let the i-th symbol be 
"1". Then 



-^m,_i-m,((-"^l - 7712 - ... - mi_i)K) = Jo(0). (60) 



The previous term is 



Jmi.2-mi_Ai-'^l -^2- ... -mi_2)K). (61) 

For both to be Jo(0) it is required that mj_i = 0, and mj_2 — mj_i = 0, implying mj_2 = resulting 
in 

^m,_3((-"^i -m2- ... - mi^3)K) = 0, 

for rrij-s 7^ 0. Therefore if the term before the i-th one is "1" ( and we have two neighbors that 
are "l"s ) either mi_3, and all rrij with j < i — 2, vanish and all factors before the i-th are "l"s, in 
contradiction with the fact that for m 7^ the sequence starts with an x, or the contribution of the 
sequence vanishes ( when one of the rrij does not vanish ) . Now one has to find the strings (^) with 
the maximal number of ones subject to given values of m, m' and q. For this purpose strings with 
alternating "x" and "1" are constructed. 

The "x" represent factors Jm,-mi+-i{—MiK) where Mi = m + I]i=i ^i and we have to choose the 
rrii so that the Jmi-mi+i{—MiK) are of maximal magnitude. Consider the string: 

Jn.,^rnU-MlK) J„,^,_^,^, (-M^+iiT) J^,^,_„,^3 (-M,+2i^) J^,+3-„.,+, (-M^+gi^) (62) 



where the factors exp(— ^m^) were omitted for the sake of brevity. Requiring that the second and 
fourth factors are "1" yields mi^i = m/+2 and m/+3 = 771^+4 as well as M/+i = M/+3 = 0. Therefore 
Ml = -m/+i and M1+2 = 'm/+2 = -mi+3. This implies m/+i = m/+2 = -m/+3 = -771^+4 = m*, and 
this string takes the form: 

Jm,-m*{m*K) Jo(0) J2m*{-m*K) Jo(0).... (63) 



Continuation of the string to the left requires mi = —m*. The factors "x" are J2m*{—fT^*K) = 
J-2m*{fn*K) = J2m*{f^*K). For cach value of K we choose the value of m* so that | J2m*{f^*K) \ is 
maximal, namely 

I J2m*('m*K) |=max| J2m{niK) I . (64) 

m, 
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Now it is left to match this string to the ends that are determined by m, m' and q. The term 
is the sum of terms of the form: 

C^^\m,m*) (^J2m4m*K)exp{-a^m*^)y' C'^'\m*,m',q) (65) 

where n' is an integer of the order n/2. The string (^) is of period 4 therefore also the end terms 
are of period 4 in n. One can always find enough values of rrii in the beginning and in the end of 
strings in ( p7D so that they take the form (|65|). The end terms C'^^\m,m*) and C^^\m*,m',q) are 
sums of the contributions of these rrii. Some of the contributions to the end terms are presented in 
App. B. The end terms do not affect the resonance. Therefore the largest resonance associated with 
the fast decaying modes, corresponding to the slowest one, is up to the fourth root of the identity: 



A J2m* {m*K) exp I — — m*2 j , (66) 

independent of m, m' and q. The reason that z is determined only up to the fourth root of the 
identity is the period 4 of the string (0). The resonances associated with the other fast decaying 
modes cannot be determined in the framework of the perturbative expansion of the present work. 

The Frobenius-Perron operator is the evolution operator U in the limit of vanishing noise. There- 
fore the Ruelle resonances are the poles of matrix elements of the resolvent R in this limit. They 
form several groups. There is 

zo = 1, (67) 

that is related to the equilibrium state found for m = m' = q = 0. The resonances corresponding to 
the relaxation modes related to the diffusion in the angular momentum space are given by 

^'^ = """"P ( ~^^^^ ~ 2J2(ir)) j . (68) 

Finally, the largest resonance related to relaxation in the 6 direction is, up to the fourth root of the 
identity. 



p2,n*{m*K), (69) 



where m* is chosen so that the expression is maximal for a given value of K. The corresponding 
relaxation rates 7^ and 7 are defined by 



Zk = e-^^ (70) 



and by 



\z\ = e-\ (71) 

leading to 7^ = \lnzk\ and 7 = |/n|5||. The last resonance may take the four values 

z = ±e-\ (72) 

and 

z = ±ie-^. (73) 

The perturbative calculation enables to compute only \z\. 
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IV. NUMERICAL EXPLORATION OF RELAXATION 

In Sec. 3 the Ruelle resonances were calculated for large K and extrapolated from finite to 
vanishing variance of the noise a. Finite noise has the effect of truncation of the matrix of the FP 
operator and the limit a — *> is the infinite matrix limit. In the limit K ^ oo complete stochasticity 
takes place, while for finite K the system is a mixed one, but for large K the chaotic component 
covers nearly all of phase space. The results of Sec. 3 were obtained as the leading terms in an 
expansion in powers of 1/yK. In the present section the results will be tested numerically for finite 
K and a = 0. The phase space (^ with various values of s will be used. The resonances of the 
type (^), corresponding to diffusion in angular momentum J and of the type ( |7TD corresponding 
to relaxation in the 9 direction will be calculated numerically from the relaxation rates of various 
perturbations to the uniform invariant density. For large s, the relaxation of the diffusion modes 
(E9) (with small k) is slow and these dominate the long time behavior. To see the angular relaxation 
modes (|7TD one has to eliminate the slow relaxation. This can be done either by the choice of small s 
or by the use of distributions that are uniform in the momentum J . Evolving an initial distribution 
g for n time steps and projecting it on a distribution / defines the correlation function: 

CfM = {f\lf-\g). (74) 

For a chaotic system, for large n it is expected to decay as 

Cf.in) ~ e-^" (75) 

and the relaxation rate 7 is computed numerically from plots of Cfg{n) as a function of n. In what 
follows the distributions g and / will be selected from the Fourier basis (|^ so that 7 is expected to 
take the values 7^ or 7. Relaxation of the form ([75| ) is expected to hold in the chaotic component. 
An efficient way to calculate correlation functions like ( [7^ ) projected on this component is from a 
trajectory in phase space. By ergodicity it samples all phase space in this component. The phase 
space integrals involved in the calculation of ( |74D are replaced by time averages along the trajectory. 
The trajectories were started in the vicinity of the hyperbolic point (vr, 0) and iterated for a large 
number of time steps, A^. It was verified for several cases that the results "equilibrize" , namely they 
do not depend on A^ for large A^. The correlation function is calculated from the formula: 

1 ^ 
Cf,{n)=\un^-Y.fiM^ + n) (76) 

1=1 

where /(j) and g{j) are the values of / and g at the j-th time step. We first calculate numerically 
the slow relaxation rates •jk ( |TDD related to diffusion and then turn to calculate and 7 ( [7I| ) related 
to relaxation in the 6 direction. 

A. The Diffusive Modes 

The relaxation rates expected from the approximate calculations of Sec. 3 for the diffusive modes 
are given by (^8|) or 

7fc = ^DiK), (77) 

where D{K) is the diffusion coefficient (^ for cr = 0. To test this relation, the correlation function 
(0) was calculated for various distributions g and / from the Fourier basis (|^) and plots like the 
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ones presented in Fig. 1 were prepared. The slope is 7^ and the values of D{K) are extracted with 
the help of (|7^) for various values of k and s. In Fig. 2 these values of D{K) are depicted for large 
values of the stochasticity parameter K. Excellent agreement with the theory is found: (a) The 
value of D is found to be independent of k and s; (b) It agrees with the theoretical prediction (|^). 
We find indeed that for long time the behavior of distributions is indeed as for a diffusive process. In 
the past it was checked that only the second moment of the momentum grows linearly as expected 
for diffusion. The effect of sticking to the accelerator modes was not observed for the values of K 
used for Fig. 2 since the size of the accelerated region is small and therefore special effort is required 
to observe these effects in numerical calculations [|I^]. These are expected to be important for 
relatively small values of K where the accelerated regions are larger. 

In Fig. 3 the correlation function is plotted for relatively small values of the stochasticity parameter 
K where larger deviations from the theory presented in Sec. 3 are expected. The diffusion coefficient 
as a function of K is presented in Fig. 4. Deviations of the numerical results from the analytical 
predictions are found for some values of K. Also for these the decay of correlations is found to be 
exponential and the diffusion coefficient extracted for all modes by (^7D is the same. Therefore the 
behavior that is found is indeed diffusive, but the value of diffusion coefficient for some values of K 
is larger than the one that is theoretically predicted. This is a result of sticking ( for finite times ) to 
accelerator modes. For most values of K the value of D found from (|77| ) agrees with the one found 
from direct evaluation of trajectories in the chaotic component. The theoretical errors (marked by 
dashed line in Fig. 4 ) were estimated from the next term of the formula of Rechester and White 
for the diffusion coefficient 
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The actual errors are larger due to nonperturbative nature of the accelerator modes and the 
surrounding regions (such modes cannot be found in an expansion in 1/Vk ) . Since in all calculations 
only trajectories belonging to the chaotic component were propagated real acceleration is avoided. 
The trajectories used in the calculation of the correlation function by ( |76D effectively generate a 
projection on the chaotic component of phase space. 
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FIG. 1. The function Cfg{n) (semilogarithmic plot) for: (a) f = g = (pio,K = 20, s = 370, ; (b) / = g = 02o, K = 30,s = 900; 
(c) / = <; = <f)50,K = 40, s — 3200; (d) / — (j>i2,g = (kn, K — 27, s = 450. The dashed hne represents the best fit to the data. 
The number of iterations is A?^ = 8 * 10®. 
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FIG. 2. The diffusion coefficient D for iT > 10 as found from plots like the ones presented in Fig.l for: the first mode (la) 
(squares), the second mode (lb) (stars), the fifth mode (Ic) (circles), correlation function for the off diagonal first mode (Id) 
(pentagram) and other off diagonal correlation functions (triangles), compared to the theoretical value (solid line). 
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FIG. 3. The function Cfg{n) for (a) / = g = 0io,iC = 7, s = 250, (b) / = g = (l>2o,K = 8, s = 510, (c) 
f ~ g ~ (psOjK — 3,s — 340, (d) / = <j}ii,g — (j>i2,K — 17, s — 225. The dashed line represents the best fit to the 
data. The number of iterations is TV = 8 * 10^. 
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FIG. 4. The diffusion coefficient D for Ti" < 20 as found from plots like the ones presented in Fig. 3 for: the first mode (la) 
(squares) , the second mode (lb) (stars), the fifth mode (Ic) (circles) and off diagonal first mode (Id) (pentagram), compared 
to the theoretical value (solid line). The dashed line represents the approximate error. The values of D obtained by direct 
simulation of propagation of trajectories are marked by diamonds. 



B. Angular Relaxation 



In order to observe the angular relaxation mode it is required that no relaxation in the J direction 
is present, because such a relaxation if present, is expected to dominate the long time limit. Since 
the results are independent on s, we use s = 1. For this purpose we take g = (pkm' (see (^) so 
that q = k/s is an integer, and / = 0om- From ( |55D and (^) one concludes that the slowest of the 
angular relaxation rates is 



7 = --ln[max(|J2m*("^*^)|)]- 



(7J 



The absolute value of the correlation function Cfg{n) is presented in Figs. 5 and 6 for (? = {/)o2 and 
/ = (pQi and for g = f = 0oi respectively for several values of K. The numerical calculations are 
complicated since the relaxation is fast, with a characteristic time of the order of one time step. 
Moreover there are oscillations of the correlation function while ( [78D is just the envelope. In Figs. 
5 and 6 the best fit to the envelope is marked by a dashed line. The slope of the dashed line is 
the numerical estimate for the relaxation rate. In Fig. 7 the numerical estimate is compared with 
the theoretical prediction. The error in the theoretical prediction is estimated as the value of the 
next order contribution to a„. This results from a term where the "l"s in sequences corresponding 
to (|63D and (|65|) is replaced by an "x" that represents a Bessel function of order 1/a/a , leading 
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to an error of the order In (1 ± IjyK) in the relaxation rate. It is difficult to estimate the error 
resulting from the numerical procedure of calculating the relaxation rates. The reason is that near 
the origin of the correlation function a large number of modes contributes. On the other hand, in 
the tail of the correlation function, where only one relaxation rate is dominant, the signal is too 
small. Nevertheless, the comparison between our numerical and theoretical results shows a good 
qualitative agreement. 
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FIG. 5. The absolute value of the function Cfg{n) for / = <^oi,P = 0O2 and (a) K=16.3, (b) K=19.5, (c) K=12, (d) K=16. 
The dashed line represents the best fit to the data. The values s = 1 and A*' — 10* were used. 
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FIG. 6. The absolute value of the function Cfg{n) for / = 3 = (fxn and (a) K=10.7, (b) K=14.3, (c) K=12.55, (d) K=14.7. 
The dashed line represents the best fit to the data. The values s = 1 and A*' — 10* were used. 
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FIG. 7. The fast relaxation rate 7 as found from plots like Figs. 5 and 6 for / = ^oi, 9 ~ <?f>02 (triangles) , / — 0o2, 
g — </f>04 (circles), f = g = (^02 (stars) and f = g = </f>oi (squares), compared to the theoretical value ( |78| ) (solid line). The dashed 
lines denote the theoretically estimated error. Here we used s — 1 and A'' = 10 . 

V. SUMMARY AND DISCUSSION 

Relaxation to equilibrium was studied for the kicked rotor that is a standard system for the 
exploration of classical chaos in driven systems and its quantum mechanical suppression. Relaxation 
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and diffusion are important concepts in statistical mechanics. Here they were studied for a mixed 
chaotic system. Very httle is known rigorously about such systems although most models describing 
real physical systems are mixed, namely in some regions of phase space the motion is regular while 
in some regions it is chaotic. 

In this work the kicked rotor was studied in a phase space that is the torus defined by (^. The 
relaxation of distributions in phase space takes place in stages. First, the inhomogeneity in 9 decays 
with rapid relaxation rates, the slowest of them is 7. Then relaxation of the inhomogeneities in the J 
direction takes place with the relaxation rates related to the diffusion coefficient via (^). Diffusion 
was previously believed to be a good approximation for the kicked rotor, but here, to our best 
knowledge, the various time scales were analyzed carefully for the first time. In particular we have 
found the time scale, I/7, below which the diffusion approximation does not hold since relaxation 
of correlations in the angle direction still takes place. 

There is a clear relation between the relaxation of inhomogeneities in 9 and the diffusion constant 
since 

n 

< {Jn+1 - Jof >= ^ K^ < sin9isin9j >, (79) 

i,j=0 

where Jj and 9i are the momentum and angle before the i-th kick. For a chaotic trajectory 

< sin9isin9j >=< sin9osin9\i^j\ >= Cff{\i — j|), (80) 

where Cff{\i — j\) is the correlation function (^) with / = 0oi- If the sum YnLoCffii) converges, 
as is the case where Cff falls off exponentially, diffusion is found and the value of the diffusion 
coefficient is 

T^2 00 

^ = ^EC//(0. (81) 



In App. C we show that ([50|) that was obtained by Rechester and White in |14| is just 



D = ^j: Cfj{i). (82) 

^ i=-2 



A derivation that is very similar is presented in [Q. If the sum diverges one obtains anomalous 
diffusion. 

Finite noise leads to the effective truncation of the evolution operator (^) . In the basis (|^ it 
means that it results in limited resolution. Moreover for a > the operator U is nonunitary. The 



approximate eigenvalues of U given by (0), that were found in this work are 1 and Zk of (^Sj) ( if 
(PH) satisfied ) and z of (|7T[). In our approximation method we could not obtain many eigenvalues 
related to angular relaxation modes. Because of the effective truncation, ip-^, the eigenfunction of f/, 
can be expanded in terms of the basis states (0). The relaxation rates of these states are —ln{zk) 
and —ln{\z\), where Zk and z are given by (^) and (|7ID. In the limit a — ;> the evolution operator 
is unitary, ip^ approach some generalized functions while Zk and z approach the values of the poles 
of the matrix elements of the resolvent R of (^5]) obtained from the extrapolation from \z\ > 1 
(corresponding to the < e ^ 0, used in the standard definition of the Green's function ). 

These are the Ruelle resonances that are related to the relaxation rates via ( ffOD and (ffTI). This is 
very similar to the situation for hyperbolic systems such as the baker map. For hyperbolic systems 
the Ruelle resonances ( related to the relaxation rates ) approach fixed values inside the unit circle 
in the complex z plane in the limit of an infinite matrix for the evolution operator or of infinitely 
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fine phase space resolution. This was found to be correct also here when one takes the limit a ^ 
in (|49D and (|66|) resulting in (^) and (|69D. Numerical tests in absence of noise confirm that the 
analytical results provide a good approximation for the relaxation to equilibrium and diffusion in 
the chaotic component. Results of similar nature were found in the standard map with s = 1 for 
some values of K ||T^, for the "perturbed cat" map [0, and also for the kicked top [^. In all 



these works it was found, within the approximations used, that the leading resonances are either 
real or form the quartet (iA, ±iA) where A is a real number satisfying < A < 1. The generality 
of this form should be subject to further research. For the kicked top it was attributed [^ to the 
dominance of an orbit of period 4. 

In mixed systems, such as the kicked rotor, even in the chaotic components there is sticking to 
regular islands and acceleration modes. Noise eliminates this sticking . The analytic formulas ( |68|) 
and (^) are obtained from an expansion in powers of 1/yK for finite variance of noise a^ and 
the limit a^ ^ is taken in the end of the calculation. A nonvanishing value of a^ assures the 
convergence of the series (|57D . Appearance of the islands and the sticking is a nonperturbative effect 
and therefore it is not reproduced in our theory. For this reason in absence of noise the results are 
only approximate. The effect of the sticking is extremely small for most values of the stochasticity 
parameter K, as verified by the numerical calculations without noise. 

The physical reason for the decay of correlations is, that in a chaotic system, because of the 
stretching and folding mechanisms there is persistent flow in the direction of functions with finer 
details namely larger \k\ and |m| in our case. Consequently the projection on a given function, for 
example one of the basis functions (|^) in our case, decays |^,29|. The crucial point is that this 



function should be sufficiently smooth. This argument should hold also for the chaotic component 
of mixed systems. In the present paper the actual relaxation rates were calculated. Here noise was 
used in order to make the analytical calculations possible. In real experiments some level of noise is 
present, therefore the results in presence of noise are of experimental relevance. It was shown with the 
help of the Cauchy-Hadamard theorem |2^ ( see discussion following ( ^2] ) ) that for s >> K >> 1 



exponential relaxation to the invariant density takes place with the rate 71 = D{K)/s'^, where D{K) 
is given by (^ . It was deduced from the radius of convergence for the series of the matrix element 
of R' (see (p9|)). This rate is independent of a. It is found for all functions that can be expanded 
in the basis (|^, with an absolutely convergent expansion. It excludes for example functions of the 
form (|22|). We believe this statement can be made rigorous by experts. 

For the baker map it was found that the resolvent of the evolution operator of the Quantum 
Wigner function, when coarse grained has the same poles as the classical Frobenius-Perron operator 



21| . We believe it should hold also here. The fact that the Ruelle resonances of the modes of slow 
relaxation are Zk that are identical to the ones of the diffusion operator gives additional support to 
approximations made for the calculation of the ensemble averaged localization length in [^. 

Finally, the Ruelle resonances, that were introduced and established rigorously for hyperbolic 
systems can be used to describe relaxation and transport in the chaotic component of mixed systems. 
Here it was demonstrated for the kicked-rotor. 
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APPENDIX A: THE MATRIX ELEMENTS OF THE EVOLUTION OPERATOR IN PRESENCE OF 

NOISE. 

In this appendix the matrix elements in the representation (0) are calculated. For this purpose 
(pop is transformed to the Fourier representation by 

(fc2m2|f/|A;imi) = dO dJ dO' dJ'{k2m2\je){je\Unoise\J'd'){J'd'\UK\kimi). (Al) 

Jo Jo Jo Jo 

Substitution of (|18D and ([T9|) yields 

(k2m2\U\kimi) = / d9 dJ dO' exp(— i??i26')exp i (A2) 

Jo Jo Jo \/27r \/2'KS \ s 



27rV27fs 

V" — exp{im{6' — 6)) exp —imJ m^ — ^^^ exp(zmi6'^) exp i — (J + KsinO') \ . 

m 27r y 2 y v27r v 27rs \ s / 

Terms containing 9 are exp(i?7i(— ^)) exp (— 277126'). Integration over 9 yields 5m.-m2 leading to 

1-27X3 r-2-K 1 / _i. /\ 

(k2m2\U\hmi) = dJ d9'^^exp{i — —] (A3) 

Jo Jo v27rs \ s J 

— exp(i(— 77126'') exp im2J 7712 ) —^= exp{imi9') exp i — (J + Ksin9') ) . 

27r \ 2 J v27rs \ s ) 

Integration over J results in 6k2-ki,m2s, yielding 

{k2m2\U\kimi) = ^ / d9' exp{-i{m2 - 7711)6'') exp I — — 7712 1 exp ii—Ksin9' j 5k2-k^,m2s (A4) 
Finally with the help of the integral representation for Bessel functions: 

Jmiz) = 7^ / (i6'exp(— i77i6') exp(i2;sm6') 

277 Jo 

one obtains (^T]) . 

APPENDIX B: THE END TERMS IN STRINGS OF THE FAST MODES 



In this appendix, possible examples for contributions to the end terms in ( [55| ) are presented. The 
left end term is a sum of terms of the form 



2 

J2m+m* {-mK) J-m-2m* {m* K) Jo(0) exp — (771^ + (771 + 771*)^ + 771*^) 

+ ^ Jm-mii-mK) J2rm+m+m*{{-m - mi)K) J^m-2fn*~mi{m*K) Jg(0) 
mi 

2 

exp im + nil + [m + m + '^ij +^12 



(Bl) 
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where in the first term mi = {m + m*) and m2 = m^ = m* while in the second term 7712 = 
— (rra + m* + m-i) and m^ = m^ = m*. For g 7^ the right end term is a sum of terms of the form 

C^'\m*,m',q) = 



^2 



exp 



'— yni*'^ + (m* + m„_i + g)^ + +m'^_-^ 



(B2) 



where mn-2 = —{m* + rrin-i + Q') and m„_3 = m*. For g = one has to take m„_i = m' and the 
end term consists of a sum over 'mn-2- 

APPENDIX C: THE RELATION BETWEEN THE DIFFUSION COEFFICIENT AND CORRELATION 

FUNCTION. 

In this appendix the relation between (|8^) and (|50D will be derived ( for a somewhat similar 
derivation see |^ ). For this purpose we note that 



/-27r rjg /■2-Ks fj J ^ 1 

^"^""^^l 2^1 ^,''^(^U-'zne = --m-l\-{OM]U-[\0,l)-\0,-l)], (CI) 

where the representation \k,m) (see (|^)) is used. The matrix elements of U are given by (^) and 
the matrix elements of If^ required for the present calculation are 

(0,m2|f/2|0,mi) = J2„, (-m2K)e-^'™'(5„„„^, (C2) 

as can be easily obtained from multiplication of two matrices of the form (0). From ( |C1| ) it is 
clear that Cff{0) = |. Inspecting ( plD with ki = k2 = one notes that it is required that also 
mi = 1712 = 0, therefore C//(l) = 0. Substitution of ( |C2| ) in ([Cl] ) yields 

Cff{2) = -^-J2{K)e-\ (C3) 



Using the fact that Cff{—n) = Cff{n), substitution of the values of C//(0) and C//(2) into ( ^2] ) 
yields the expression ( |50D that was obtained by Rechester and White [|I^. Because of the discussion 
following ( pU] ) the correlation functions Cff (n) with n > 2 lead to terms that are of higher orders in 



l/VK than §^) 
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